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Detecting  Periodic  Components  in  a  White  Gaussian  Time  Series 


By  Shean-Tsong  Chiu 
Rice  University 


SUMMARY 

A  family  of  tests  for  periodic  components  in  a  white  Gaussian  series  is  proposed.  The  test  is 
based  on  a  statistic  which  is  proportional  to  the  ratio  of  the  maximum  periodogram  to  the 
trimmed  mean  of  the  periodograms.  The  asymptotic  distribution  of  the  statistic  is  obtained.  It  is 
shown  that  the  test  proposed  and  Fisher’s  test  have  the  same  asymptotic  powers  at  the  alterna¬ 
tive  hypotheses  that  the  series  contains  a  single  periodic  component  at  a  non-zero  Fourier  fre¬ 
quency.  The  tests  are  applied  to  detect  the  eigenfrequencies  of  the  Earth.  The  proposed  test 
detects  some  peaks  which  Fisher’s  test  fails  to  detect. 

Keywords:  Detecting  periodic  components;  Fourier  transform;  Periodogram;  Eigenvibration  of  the 
Earth. 

1 .  Introduction 

In  this  study,  we  consider  the  situation  where  we  have  T  consecutive  observations 
X(f),<=0,l,  .  .  .  ,  T-l  on  a  time  series.  We  are  interested  in  detecting  periodic  components  in 
the  series.  The  periodogram  has  been  used  to  search  for  hidden  periodicities  and,  naturally,  it 
became  a  tool  for  testing  the  presence  of  periodic  components.  Fisher  (1929)  introduced  a  test  for 
periodicity,  which  is  based  on  the  periodograms  at  the  Fourier  frequencies;  the  test  statistic  is  the 
ratio  of  the  maximum  of  the  periodograms  to  the  total  sum  of  the  periodograms.  Fisher  also  gave 
the  exact  distribution  of  the  test  statistic  under  the  null  hypothesis  that  the  series  is  a  white 
Gaussian  series.  The  asymptotical  distribution  of  the  test  statistics  was  given  in  Bloomfield 
(1976).  Fisher’s  test  was  used  in  Nowroozi  (1965),  (1966)  for  detecting  the  eigenfrequencies  of  the 
Earth.  The  test  was  also  mentioned  in  many  books  concerning  time  series  analysis.  Hartley  (1949) 
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proposed  a  test  based  on  the  ratio  of  the  maximum  periodograms  to  the  residual  sum  of  squares, 
where  the  residuals  are  obtained  by  fitting  a  sinusoidal  wave  to  the  series. 

In  practice,  we  might  suspect  that  the  series  consists  of  several  periodical  components;  in 
this  case,  the  periodogram  will  have  multiple  peaks.  Multiple  peaks  could  also  be  caused  by  the 
leakages  from  nearby  peaks,  by  the  harmonics  of  a  fundamental  frequency  or  by  a  decaying 
sinusoidal  wave  such  as  the  eigenvibration  of  the  Earth  (Bolt  and  Brillinger  (1979)).  In  these 
situations,  we  are  often  not  satisfied  just  with  rejecting  the  null  hypothesis  of  no  periodic  com¬ 
ponent  but  are  also  interested  in  detecting  all  significant  peaks  in  the  periodograms.  As  indicated 
by  Whittle  (1954),  Fisher’s  test  is  quite  conservative  in  detecting  multiple  peaks.  Shimshoni 
(1971)  also  observed  that  the  method  was  unable  to  detect  some  periodic  components  in  the  study 
of  the  eigenvibrations  of  the  Earth. 

For  detecting  the  presence  of  multiple  peaks,  various  methods  were  proposed  and  studied. 
Some  of  the  methods  were  developed  for  testing  the  outliers  in  exponential  samples.  We  remark 
that  these  two  testing  problems  are  closely  related.  Whittle  (1954)  suggested  removing  the 
significant  periodic  components  from  the  series  and  applying  Fisher’s  test  to  the  residual  series.  It 
might  not  be  easy  to  do  so  when  the  periodic  components  are  not  pure  sinusoid  waves,  such  as  in 
the  case  of  sinusoid  waves  with  non-constant  amplitudes  or  phases.  Shimshoni  (1971)  suggested 
using  the  ratio  of  the  r-th  largest  periodogram  to  the  sum  of  periodograms.  The  main  difficulties 
of  this  approach  are:  (l)  This  method  implicitly  assumes  that  the  frequencies  of  the  periodic  com¬ 
ponents  are  at  the  Fourier  frequencies;  this  is  usually  not  true  in  practice.  (2)  The  number  of  the 
periodic  components  must  be  known.  (3)  If  the  number  of  the  periodical  components  is  smaller 
than  r,  then  this  test  will  have  very  small  chance  to  reject  the  null  hypothesis  even  when  the 
series  does  contain  some  significant  periodic  components.  Lewis  and  Fieller  (1979)  studied  another 
test  which  is  based  on  the  ratio  of  the  sum  of  the  r  biggest  periodograms  to  the  total  sum  of  the 
periodograms.  This  test  will  reject  the  null  hypothesis  when  there  is  only  one  huge  peak  in  the 
periodograms.  Therefore,  the  test  is  not  able  to  tell  us  which  peaks  are  significant.  Another  disad¬ 
vantage  of  this  test  is  that  it  loses  significant  amounts  of  power  when  the  number  of  periodic  com- 
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ponents  contained  in  the  series  is  different  from  r . 

Likes  (1966)  and  Kabe  (1970)  studied  the  test  statistic  which  is  the  ratio  of  the  difference 
between  the  biggest  and  the  r-th  biggest  periodograms  to  the  biggest  periodogram.  Siegel  (1980) 
proposed  a  method  based  on  the  sum  of  the  periodograms  (normalized  by  the  total  sum  of 
periodograms)  which  are  bigger  than  a  portion  of  the  critical  value  of  Fisher’s  test.  One  must 
select  a  proper  portion  in  order  to  apply  this  test.  The  best  choice  is  dependent  on  (and  is  sensi¬ 
tive  to)  the  number  of  periodic  components  in  the  series.  These  two  methods  also  share  some  of 
the  difficulties  of  the  tests  mentioned  earlier. 

In  the  next  section,  we  describe  some  basic  notations  and  properties  of  the  Fourier  transform 
and  periodograms.  Fisher’s  test  and  its  asymptotical  distribution  are  reviewed  in  Section  3.  In 
Section  4,  we  propose  a  test  statistic  which  is  proportional  to  the  ratio  of  the  biggest  periodogram 
to  the  trimmed  mean  of  the  periodograms.  The  asymptotic  distribution  of  the  test  statistic  is  also 
derived.  In  Section  5,  we  derive  the  asymptotic  power  functions  of  the  tests  against  the  alterna¬ 
tive  hypothesis  that  the  series  contains  a  sinusoidal  wave  at  one  of  the  Fourier  frequencies.  We 
show  that,  asymptotically,  the  test  proposed  has  the  same  power  as  Fisher’s  test  has.  In  Section 
6,  we  apply  the  proposed  test  to  detect  the  eigenfrequencies  of  the  Earth’s  free  vibrations.  We 
show  that  our  method  can  detect  some  of  the  peaks  which  Fisher’s  method  fails  to  detect.  The 
asymptotic  results  are  supported  by  simulations. 

In  this  research,  we  are  concerned  with  the  case  that  the  underlying  series  is  a  white  series. 
The  case  of  non-white  series  was  considered  in  Hannan(1961)  and  Whittle(1952,  1954). 

2.  Basic  notation  and  properties  of  periodogram 

The  data  we  consider  here  are  a  sequence  of  observations  taken  at  equal  time  intervals.  The 
series  consists  of  a  deterministic  series  and  a  noise  series: 

X(t)  =  ?+S(t)+c(i)  1=0,1,  .  .  .  ,  T-l..  (2.1) 

Here  e(t)  is  a  white  Gaussian  series  with  means  zero  and  variances  a 2  and  the  deterministic  part 

K 

5(t)  =  £  A*cos  [lkt+4>k) 

k  =  1 


(2.2) 


-  4  - 


consists  of  K  sinusoidal  waves  at  frequencies  with  amplitudes  A*  and  phases  <£*;  the  fre¬ 

quencies  are  not  necessarily  the  Fourier  frequencies.  Without  loss  of  generality,  we  assume  p=0. 
We  are  interested  in  testing  for  the  existence  of  S(t).  The  null  hypothesis  can  be  formulated  as 
H0  :  5(l)  =  0.  X(t)  is  a  Gaussian  white  series  under  H0 

The  Fourier  transform  of  the  series  X(t)  at  frequency  ui  is 

dJrr)M  =  X)X(<)exp(-twl)  .  (2.3) 

«= o 

The  Fourier  transforms  are  usually  evaluated  by  the  fast  Fourier  transform  at  the  Fourier  fre¬ 
quencies  u)j=2irj/T.  The  Fourier  transforms  at  the  Fourier  frequencies  are  referred  to  as  the 
discrete  Fourier  transform  of  X(f).  The  Fourier  transforms  d[U  and  d)r*  have  corresponding 
definitions.  The  periodogram  of  the  series  is  defined  by 

7x(w)  =  (2.4) 

The  periodograms  and  Is (w)  are  defined  in  a  similar  way.  Under  the  null  hypothesis, 

7x(w)  =  7'(cj)  and  the  periodograms  7x(wy )  at  the  Fourier  frequencies,  uij=2/nj/T, 
j  — 1,2,  .  .  .  ,(T/2)-l,  where  T  is  an  even  number  (for  odd  T,  j  — 1,2,  .  .  .  ,  [T/2]),  are  identi¬ 
cally  independently  distributed  according  to  the  exponential  distribution  with  mean  \=(T2/(2n). 

For  the  series  with  deterministic  part  S(t)^0,  djfT^(w)  =  df^(w)  +  d£r).  2 7a(w)/<t2  has  a 
noncentral  chi-square  distribution  with  two  degree  of  freedom;  the  noncentrality  parameter  is 
21s (lo)/ o2 .  For  detail  discussion  of  the  properties  of  the  Fourier  transform  and  periodograms,  the 
readers  are  referred  to  Brillinger  (1975). 

3.  Fisher’s  test  and  its  asymptotic  distribution 

For  testing  the  existence  of  periodic  components,  Fisher  (1929)  proposed  a  test  based  on  the 
test  statistic 

ft  =  max  Ix(wy)/  £ Ix(«y)  •  (3.1) 

1<J<»  /=! 

Here  n  —  (T /2)-l  for  even  T  and  n  =  (T- 1)/2  for  odd  T.  In  that  paper,  Fisher  gave  the  exact 
distribution  of  Fj  under  the  null  hypothesis.  Bloomfield  (1976)  derived  the  asymptotic 
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distribution  of  the  statistic  Zj  —  nFj  -  logn , 

P (ZT<z)  ~  exp{-exp(-z)}  .  (3.2) 

Since  X  =  n_I^7(wy)  converges  to  X  in  probability  and  X  has  variance  of  order  n-1, 

P(ZT<z)  =  P{Z'T<z)  +  0(n-1/2logn)>  (3.3) 

where  z  is  any  fixed  constant  and  Z'T—  max{7x(wy)/X}  -  logn.  Also 

P  (Z'T<z)  =  l-[l-exp(-2)/n]*  =  exp{-exp(-z)}  +  0(n_1).  (3.4) 

So  we  can  write  (3.2)  more  precisely, 

P(ZT<z)  =  exp{— exp(— 2)}  +  0(n~1//2logn).  (3.5) 

Let  I{uij)=Ix{u,-)/  X  be  the  normalized  periodograms,  then  nFT  =  max7(u>y).  We  see  that 
Fisher’s  test  statistic  is  proportional  to  the  maximum  of  the  normalized  periodograms.  When  the 
series  contains  periodic  components,  the  periodograms  at  frequencies  close  to  the  frequencies  of 
these  components  have  large  magnitude  and  these  periodograms  can  be  viewed  as  outliers  in  an 
exponential  sample.  From  this,  it  is  easy  to  understand  the  weakness  of  Fisher’s  test.  Fisher’s  test 
uses  the  sample  mean  to  normalize  the  periodograms;  however,  the  sample  mean  is  sensitive  to 
outliers.  The  sample  mean  tends  to  be  larger  than  the  true  mean  when  the  series  contains  periodic 
components;  therefore,  7(wy)  tends  to  be  smaller  than  7(w;)/X.  So,  the  power  of  Fisher’s  test  will 
be  smaller  than  the  power  of  the  test  based  on  the  maximum  of  7(u>y)/X.  This  suggests  that  we 
consider  tests  based  on  the  maximum  of  periodograms  normalized  by  a  robust  estimate  of  X. 

4.  The  proposed  test  and  its  asymptotical  distribution 

The  test  proposed  here  based  on  the  test  statistic 

M 

RtW  =  iJY,U.  (4.1) 

»  =  1 

where  n  is  defined  in  Section  3,  and  7^72<  •  •  •  <7,  are  the  order  statistics  of  the  periodo¬ 
grams.  f3  is  a  constant  between  0  and  1  which  determines  the  proportion  of  periodograms 
trimmed.  In  most  case,  /?=0.9  or  0.95  will  give  good  result.  The  performance  of  the  test  is  not 
very  sensitive  to  the  choice  of  (3.  Letting  /?=!,  we  obtain  Fisher’s  test  as  an  extreme  case.  In  the 
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following  discussion,  we  assume  0<^<1  and  use  the  notation  Rj  instead  of  The  denomi¬ 

nator  in  (4.1)  is  proportional  to  the  trimmed  mean 


*,=  [«r^-  • 

*■*1 

It  is  well  known  (cf  Lehmann  (1983))  that  converges  in  mean  square  to 


(4.2) 


Jo  (  )Yexp(_a:/X^dz  =  cX>  (4-3) 

where  £(/9)  =  -log(l-/S)  and  c  =  fl+(l-0)log(l-/3).  Therefore,  Xff/e  is  an  asymptotically  unbiased 
estimate  of  X.  The  exact  unbiased  factor  can  be  found  in  Kimber  (1983). 

Following  the  argument  of  the  last  section  and  letting  ZT  =  c[n/3]RT  -  logn,  we  have 


P  (ZT<z)  =  exp{-exp(-z)}  +  O{{n0)  1//2logn)  (4.4) 

asymptotically.  This  shows  that  Zj  and  Zt  have  the  same  asymptotic  distribution  in  the  null 
case. 


Similarly,  we  can  obtain  the  asypmtotic  distribution  of 

f/r(r)=  /,,_r+i/ (4.5) 

i=i 

which  is 


P{nI/j-(r}-log(n-r  +  l)<  z}  =  exp{-exp(-z)}  +  0  {n_1/,2log(n-r  +  l)}.  (4.6) 

The  test  Uf(r)  was  used  for  testing  the  existence  of  r  periodic  components  Shimshoni  (1971). 
The  exact  null  distribution  of  Uj{r)  was  derived  in  Grenander  and  Rosenblatt  (1957).  However, 
from  the  argument  in  Section  3,  it  is  clear  that  we  should  at  least  exclude  the  r-1  biggest 
periodograms  from  the  denominator  of  (4.5).  This  leads  us  to  consider 

Vr(r)=/.-r+1/’E1/o  (4-7) 

i  =  i 


which  has  asymptotic  distribution 

P{nVrr(r)-log(n-r  +  l)<z}  =  exp{-exp(-z)}  +  O  {(n-r  +  l)-1/2log(n-r  +  l)}. 


(4.8) 
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5.  Asymptotic  power  functions  of  the  tests 

In  this  section  we  derive  the  asymptotic  power  functions  of  the  tests  described  in  the  previ¬ 
ous  sections.  Without  lost  of  generality,  we  assume  X=1  in  the  following  discussion.  We  first 
consider  the  alternatives 

Ht:  S(t)=acos(~i0t).  (5.1) 

Here  70=2;r;0/7’,  l<.;0<n  is  an  arbitrary  integer.  The  discrete  Fourier  transform  of  S(t)  is 
equal  to  aT/2  at  frequency  70  and  is  equal  to  0  at  other  frequencies.  This  implies  that  the 
periodograms  at  frequencies  other  than  70  have  the  same  distribution  under  both  hypotheses. 
Under  the  alternative  hypothesis,  Ha,  2IX(^0)  has  a  noncentral  chi-square  distribution  with  2 
degrees  of  freedom  and  the  noncentrality  parameter  2A=2a2T /Sn.  The  asymptotic  power  func¬ 
tions  of  the  tests  are  obtained  in  the  following  theorem.  Because  /„  increases  at  rate  logn,  we  let 
A  increase  at  rate  logn  . 

Theorem.  The  power  function  of  the  test  Rr  and  FT  against  the  alternatives  H,  with 
A  =J3+logn  =  c2T /87t  are,  asymptotically, 

o+(l-a)P{/x(-70)>n/0},  (5.2) 

where  a  is  the  level  of  significance  of  the  tests,  n/0=-log(-log(l-a))+logn  and  2 /x('y0)  has  a  non- 
central  chi-square  distribution  with  two  degrees  of  freedom  and  noncentral  parameter  2A . 

Proof.  Let  /0  and  r0  be  the  respective  critical  values  of  the  test  FT  and  RT  at  significance  level 
a,  then  n/0  =  ft/0+  O(n~'/2logn)  and  r0c[n/?]  =  r0e  [n/?]+ 0  {(n/?)_1/2logn  },  where 

**/o  —  «/o  +  O{n^l2\ogn)  =  r0e  [n£]  =  -log{-log(l-a)}.  (5.3) 

Since 

P{/A(7o)<*}  >  P{/t(7o)<P+logn-s}  =  1  +  0(l/n),  (5.4) 

for  any  fixed  constant  s ,  the  probability  that  /( 70)  has  a  rank  (among  all  periodograms)  smaller 
then  [n/?]  is  of  order  0(l/n)  and 

M  M 

Pffi/f=S//)  =  i+o(i/n). 

•=i  f=i 


(5.5) 
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From  the  same  argument  in  deriving  the  asymptotic  null  distributions,  we  have 

P(Rt  >  >o)  =  P(/*  >  roe[n0\  =  n/0)  +  O(n‘1/2logn).  (5.6) 

Since  P{/.‘=/'(Tf0)}=l/n  is  negligible,  (5.6)  is  equal  to 

P{K>nfo\I&I‘(lo)}  +  P{/XHo)>»»/o  and  I?<nf  0\  /„V/'bo)}  +  0(n-1/2logn).  (5.7) 

Conditional  on  I&1\' To),  the  distribution  of  /«  is  equal  to  the  distribution  of  the  maximum  of 
n-1  variables.  So,  the  first  conditional  probability  in  (5.7)  is  equal  to  a+0(l/n).  The  second 
conditional  probability  is  equal  to 

C°!oPVX{lo)>nfo  1  I‘^o)=r}fl(r\z)fQ{z)drdz,  (5.8) 

where  f  i{r  \  z)  is  the  conditional  density  of  Ic{lo)  given  /»=z  and  (70) 7^-4' >  ant^  / si2)  is  the 
density  of  I£.  Now  (5.8)  is  equal  to 

Jjog/B70IOgBp{/X('7°)>n/o  I  I<(lo)=r}f1(r\z)f^z)drdz  +  0(l/n),  (5.9) 

=  f  f°J  *  P{IX(lo)>nf0  |  Ic{lo)=r}f  i(r)f2(z)drdz  +  0(n'1/2logn), 
where  /i(r)  is  the  unconditional  density  of  Jt('y0)-  The  last  equation  holds  because 
/j(r  |  z)  —  fi(r)  4-  0(l/n)  when  z>logn.  Therefore,  the  power  function  of  the  test  RT  is  equal 
to 

a  +  (l-a)P{/x(qf0)>n/o}  +  O(n~1/2logn).  (5.10) 

The  power  of  Fisher’s  test  is 

P(Fr>/o)  =  P{/x>(n-l)/0+/o^bo)}  +  O^-^logn).  (5.11) 

Since  E[/x(70)]2=O  {(logn)2},  (5.11)  is  equal  to 

P(Fr>/0)  =  P(/,X>n/0)  +  0{n^\ogn).  (5.12) 

Therefore,  the  test  RT  and  Fisher’s  test  have  the  same  asymptotic  power  indicated  in  the 
theorem. 


We  now  consider  the  problem  of  detecting  r  periodic  components.  We  define  the  critical 
bound  of  a  test  ( FT ,  RT  or  VT)  at  significance  level  a  as  the  product  of  the  denominator  of  the 
test  and  the  critical  value  of  the  test  at  significance  value  a.  We  say  that  a  test  detects  a  peak  at 
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significance  level  a  if  the  periodogr&m  at  the  peak  is  bigger  than  the  critical  bound  of  the  test  at 
significance  level  a.  According  to  these  definitions,  FT  (or  RT)  will  reject  the  null  hypothesis  if 
and  only  if  the  test  detect  any  peak  in  the  periodogram.  The  test  VT(r)  will  reject  the  null 
hypothesis  if  and  only  if  VT(r)  can  detect  at  least  r  peaks.  The  alternatives  we  now  consider  are 

H ■  S{t)  =  J]a.cos(-7,<),  (5.13) 

1—1 

where  r  is  a  fixed  number  and  *  =1,  .  .  .  ,r  are  distinct  non-zero  Fourier  frequencies. 
Without  loss  of  generality,  we  can  assume  oj>  •  •  •  >or.  Similar  to  the  discussion  in  the  proof, 
it  is  easy  to  see  that  the  probability  of  detecting  the  r-th  peak  is  a+(l-a)P{/x('7r)>n/0},  which 
is  not  affected  by  the  r-1  biggest  peaks.  In  order  to  compare  the  probability  with  the  power  of 
the  test  VT{r),  we  assume  au  ...  ,  ar_!  are  big  enough,  such  that  the  probability  of  missing  the 
r-1  peaks  is  negligible.  Supposing  a?T /8tt  =  Br+\ogn,  we  can  show  that  the  power  of  the  test 
VT(r)  is  also  asymptotically  equal  to  o+(l-a)P{/x (^r)> n/0}. 

Remark  1.  For  the  sake  of  comparison,  we  only  consider  the  alternatives  of  pure  sinusoidal  waves 
at  Fourier  frequencies.  However,  as  shown  in  the  next  section,  the  test  RT  is  much  better  than 
the  tests  FT  and  VT  in  situations  where  the  periodic  components  are  not  at  Fourier  frequencies 
and  that  the  amplitudes  or  phases  are  not  constant. 

Remark  2.  Fr(r)  might  accept  the  null  hypothesis  while  there  are  some  (<r)  peaks  bigger  than 
the  critical  values.  This  will  cause  a  lot  of  confusion  in  interpreting  the  result. 

8.  Detecting  the  Eigenvibrations  of  the  Earth 

The  example  we  used  here  is  the  time  series  obtained  by  the  long-period  pendulums  at  the 
Grotta  Gigante,  Trieste,  following  the  1960  Chile  earthquake  (Bolt  &  Marussi  (1962)).  The  sam¬ 
pling  interval  is  2  minutes  and  the  length  of  the  series  is  2548.  Bolt  and  Brillinger  (1979)  modeled 
the  periodic  components  as  exponentially  decaying  sinusoidals.  16  sets  of  estimates  of  the  eigenfre- 
quencies,  damping  constants,  amplitudes  and  phases  were  obtained.  The  first  1200  values  of  the 
series  are  plotted  in  Figure  1.  It  shows  that  most  of  the  vibrations  die  out  after  20  hours.  We 
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compute  the  periodograms  of  the  first  600  observations.  The  critical  bounds  at  significance  level 
0.01  of  the  tests  F Rj-  and  V ?  are  also  computed.  The  results  are  shown  in  Figure  2.  The 
heights  of  the  four  horizontal  lines  in  the  plot  correspond  (from  top  to  bottom)  to  the  critical 
bounds  of  the  tests  FT,  VT(S),  i?r(0.8)  and  RT( 0.9).  We  see  that  the  tests  RT( 0.8)  and  /?r(0.9) 
both  detect  8  peaks,  while  Fj  only  detects  4  peaks.  The  8  peaks  detected  by  R?  were  also 
confirmed  in  Bolt  and  Brillinger  (1979).  There  are  only  6  peaks  bigger  than  the  critical  value  of 
Vr(8),  so  the  null  hypothesis  of  no  periodic  components  is  accepted  at  significant  level  0.01,  and 
we  are  not  able  to  make  inference  based  on  Fj-(8)  about  the  significances  of  the  peaks  which  are 
bigger  than  the  critical  bound  of  ^(8).  The  nice  feature  we  mentioned  earlier  that  the  test  R T 
is  not  sensitive  to  the  value  0  used  is  also  demonstrated  in  Figure  2.  We  see  that  the  critical 
bounds  of  /?r(0.8)  and  i?j-(0.9)  are  very  close. 

We  also  compute  the  periodograms  of  the  first  1024  observations  and  plot  them  in  Figure  3. 
The  periodograms  at  very  low  frequencies  become  dominant;  this  might  be  expected  from  Figure 
1.  The  low  frequency  components  are  primarily  caused  by  tidal  effects.  The  top  and  bottom  lines 
in  Figure  3  correspond  to  the  critical  bounds  of  FT  and  /?r(0.9),  respectively.  In  this  case,  both 
tests  detect  the  peak  at  very  low  frequency.  However,  RT  still  detects  the  same  8  peaks  detected 
earlier  while  FT  fails  to  detect  any  of  these  peaks.  This  shows  a  big  advantage  of  the  test  RT  in 
that  the  ability  of  detecting  a  particular  peak  is  not  affected  much  by  the  presence  of  other  peaks. 

Remark  3.  We  should  point  out  that  Nowroozi  (1966)  has  removed  tidal  effects  before  computing 
the  periodograms.  So,  Fisher’s  test  was  not  affected  by  the  low  frequency  component.  The  low 
frequency  components  might  not  be  easy  to  remove  in  some  cases. 

Remark  4.  The  test  Rj  failed  to  detect  some  peaks  which  were  found  in  Bolt  and  Brillinger 
(1979).  The  failure  is  caused  by  the  fact  that  periodic  components  are  decaying.  Currently,  there 
is  no  test  designed  for  detecting  decaying  sinusoidal  waves.  A  test  procedure  which  takes  the 
damping  factor  into  consideration  might  be  more  powerful  in  detecting  decaying  sinusoidal  waves. 
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7.  Discussion 

In  this  paper,  we  proposed  a  method  for  testing  the  existence  of  periodic  components  in  a 
white  Gaussian  series.  This  test  procedure  can  also  be  used  to  detect  the  outliers  in  an  exponential 
sample.  Some  of  the  advantages  of  the  proposed  test  were  discussed  in  previous  sections.  One 
small  disadvantage  is  that  we  only  obtained  an  approximation  for  the  the  exact  null  distribution. 
We  have  shown  that  the  error  caused  by  using  the  asymptotic  distribution  is  of  order  n_1,,2logn. 
The  size  of  the  series  encountered  in  practice  is  often  quite  large  and  the  approximation  becomes 
very  accurate.  Asymptotic  approach  has  an  advantage  in  that  the  asymptotic  distribution  is  easy 
to  find  and  it  is  easy  to  compute  critical  values  and  p-values. 

Acknow  ledgement 

I  would  like  to  thank  Professor  David  Brillinger  for  providing  the  data.  This  research  was 
supported  in  part  by  ARO  Grant  DAAG  29-85-K-0212  and  ONR  Grant  N  0001-485-K-0100. 

References 

[1]  Bloomfield,  P.  (1976)  Fourier  Analysis  of  Time  Series:  An  Introduction,  J.  Wiley,  New  York. 

[2]  Bolt,  B.  A.  and  Brillinger  D.  R.  (1979)  Estimation  of  uncertainties  in  eigenspectral  estimates 
from  decaying  geophysical  time  series,  Geophys.  J.  R.  astr.  Soc.,  59,  593-603. 

[3]  Brillinger,  D.  R.  (1975)  Time  Series  Data  Analysis  and  Theory,  Holt,  Rinehart,  and  Winston, 
New  York. 

[4]  Fisher,  R.  A.  (1929)  Tests  of  significance  in  harmonic  analysis,  Proc.  Roy.  Soc.  Ser.  A,  125, 
54-59. 

[5]  Grenander,  U.  and  Rosenblatt,  M.  (1957)  Statistical  analysis  of  stationary  time  series,  J. 
Wiley,  New  York. 

[6]  Hannan,  E.  J.  (1961)  Testing  for  a  jump  in  the  spectral  function,  J.  Roy.  Statist.  Soc.  Ser. 
B,  23,  394-404. 


-  12- 


[7]  Hartley,  H.  O.  (1949)  Tests  of  significance  in  harmonic  analysis,  Biomctrika,  38,  194-201. 

[8]  Kabe,  D.  G.  (1970)  Testing  outliers  from  an  exponential  population,  Metrika,  16,  15-18. 

[9]  Lehmann,  E.  L.  (1983)  Theory  of  Point  Estimation,  J.  Wiley,  New  York. 

[10]  Lewis,  T.  and  Fieller,  N.  R.  J.  (1979)  A  recursive  algorithm  for  null  distributions  for 
outliers:  I.  Gamma  samples,  Technometrics,  21,  371-376. 

[11]  Likes,  J.  (1966)  Distribution  of  Dixon’s  statistics  in  the  case  of  an  exponential  population. 
Metrika,  11,  46-54. 

[12]  Kimber,  A.  C.  (1983)  Trimming  in  gamma  samples,  Appl.  Statist.,  32,  7-14. 

[13]  Nowroozi,  A.  A.  (1965)  Eigenvibrations  of  the  Earth  after  the  Alaskan  earthquake,  J.  Geo- 
phys.  Res.,  70,  5145-5156. 

[14]  Nowroozi,  A.  A.  (1966)  Terrestrial  spectroscopy  following  the  Rat  Island  earthquake,  Bull, 
seism.  Soc.  Am.,  66, R  1269-1288. 

[15]  Siegel,  A.  F.  (1980)  Testing  for  periodicity  in  a  time  series,  J.  Amer.  Statist.  Assoc.,  75, 
345-348. 

[16]  Shimshoni,  M.  (1971)  On  Fisher’s  test  of  significance  in  harmonic  analysis,  Geophys.  J.  R. 
astr.  Soc.,  23,  373-377. 

[17]  Whittle,  P.  (1952)  The  simultaneous  estimation  of  a  time  series  harmonic  component  and 
covariance  structure,  Trahajos  de  Estadistica,  3,  43-57. 

[18]  Whittle,  P.  (1954)  The  statitical  analysis  of  a  seiche  record,  Sears  Foundation  Journal  of 


Marine  Research,  13,  76-100. 


Figure  1.  The  first  1200  observations  of  the  series,  the  sampling  interval  is  2  minutes. 


Figure  2.  The  Periodograms  of  the  first  600  observations,  the  heights  of  the  horizontal  lines 
correspond  to  the  critical  bounds  of  the  tests  Fy,  Vr(8),  Ry(0.8)  and  i?r(0.9). 

Figure  3.  The  periodograms  of  the  first  1024  observations  (in  loglO  scale),  the  horizontal  lines 
correspond  to  the  critical  bounds  of  the  tests  Fy  and  ^^(O.Q). 


t  in  Hours 


1000 


